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1. The crystal structure of Cas9/R-loop complex captured in the precleavage state (PDB ID: 5F9R) was selected as the starting model for the simulations 


performed in this study. 


2. The Mg2* ions were treated with the ion model developed by Liao el al. (DOI: 10.1021/acs.jpclett.7b02358) and the midC4 parameter set was considered 


here. The TIP4PEw water model was used to describe the water molecules, as the above Mg2* ion parameters were developed against this solvent model. 

3. Additionally, the following files (Mg_midC4.fremod and Mg.lib) are needed to be read while generating the parameter and topology files for MD simulations 
via the AmberTools tleap: 
BEM 


loadoff 


loadAmberParams Mg_minC4.fremod 


The content in Mg_minC4.frcmod: 


Mg.lib 


ES 
Mg2+ dummy model 


MASS 
MG 6.305 
Da 3.000 
Db 3.000 
Dc 3.000 
BOND 


MG-Da 800.0 0.900 
MG-Db 800.0 0.900 
MG-Dc 800.0 0.900 


ANGLE 
Da-MG-Da 
Db-MG-Db 
Dc-MG-Dc 
Da-MG-Db 
Da-MG-Dc 
Db-MG-Dc 
DIHE 
IMPROPER 
NONBON 
MG 0.706 
Da 0.000 
Db 0.000 


125.000 
125.000 
125.000 
125.000 
125.000 
125.000 


26.8375 
0.0000 
0.0000 


Dc 0.000 0.0000 
ES 


The content in Mg.lib: 
BERS ee 
'lindex array str 


"DMG" 


‘entry.DMG.unit.atoms table str name str type int typex int resx int flags intseq int elmnt dbl chg 
"MG" "MG" 0 1124972 1 19-41 ANNANN 


180.000 
180.000 

180.000 
90.000 
90.000 
90.000 


"DM1" "Da" 0 1 131072 2 1 0.500000 

"DM2" "Db" 0 1 131072 3 1 0.500000 

"DM3" "Dc" 0 1 131072 4 1 0.500000 

"DM4" "Db" 0 1 1310725 1 0.500000 

"DM5" "Dc" 0 1 131072 6 1 0.500000 

"DM6" "Da" 0 1 1310727 1 0.500000 
'entry.DMG.unit.atomspertinfo table str pname str ptype int ptypex int pelmnt dbl pchg 
"MG" "MG"0 120.0 

"DM1""Da"0 10.0 

"DM2""Db"0 10.0 

"DM3""Dc"0 10.0 

"DM4""Db"0 10.0 

"DM5""Dc"0 10.0 

"DM6""Da"0 10.0 

!entry.DMG.unit.boundbox array dbl 

-1.000000 

0.0 

0.0 

0.0 

0.0 

‘entry.DMG.unit.childsequence single int 

2 

!entry.DMG.unit.connect array int 

0 

0 

entry.DMG.unit.connectivity table int atom1x int atom2x int flags 
121 

131 

141 

151 

161 

MA 

‘entry.DMG.unit.hierarchy table str abovetype int abovex str belowtype int belowx 
"U"O"R" 1 

"R"'1"A"1 

"R"1"A"2 

"R'1"A"3 

"R"1"A"4 

"R'1"A"S 

"R'1"A"6 

"R'1"A"7 

!entry.DMG.unit.name single str 

"DMG" 

‘entry.DMG.unit.positions table dbl x dbly dblz 
0.000 0.000 0.000 

0.000 0.000 0.800 

1.000 0.000 0.000 

0.000 1.000 0.000 

-1.000 0.000 0.000 

0.000 -1.000 0.000 

0.000 0.000 -0.800 
‘entry.DMG.unit.residueconnect table intc1x intc2x int c3x int c4x int c5x int c6x 
000000 

entry.DMG.unit.residues table strname intseq int childseq int startatomx str restype int imagingx 
"DMG" 181"?"0 
!entry.DMG.unit.residuesPdbSequenceNumber array int 
0 

‘entry.DMG.unit.solventcap array dbl 

-1.000000 

0.0 

0.0 

0.0 

0.0 

lentry.DMG.unit.velocities table dbl x dbl y dblz 
0.0 0.0 0.0 

0.0 0.0 0.0 

0.0 0.0 0.0 

0.0 0.0 0.0 

0.0 0.0 0.0 

0.0 0.0 0.0 

0.0 0.0 0.0 


BEES ee 

3. And use parmed to modify the system prmtop file (named complex_wat.prmtop here, for example) derived above to include the midC4 parameter for Mg2+ 
ions: 

parmed -i parmed.in 

parmed.in contains the relevant settings read by parmed: 
BERS ES 

parm complex_wat.prmtop 

loadRestrt complex_wat.inpcrd 

setOverwrite True 

add12_6_4 @MG c4file minC4file.dat 

outparm complex_wat_Liao.prmtop complex_wat_Liao.inpcrd 
FES 


The content in minC4file.dat: 
BEM ee 
Mg-1 67.0 
BEM ee 
Note: remember to set "j1264=1" in the MD configuration file to make the Liao Mg2+ ion model effective in MD simulations! 
4. The Mg2* ions were docked at the interface between the RuvC nuclease center and the non-target DNA by referencing the crystal structure of a bacterial 
RNase H with an RNA/DNA hybrid (DOK 10.1016/.cell.2005.04.024). 
5. The built system (wild type Cas9 with a NGG PAM) underwent stages of energy minimization (2000 steps), slow heating (over 50 ps), equilibration (over 50 
ns) to reach stable. A typical configuration file for ensuing production simulation is provided below. 
SS 
&cntrl 
imin=0, 
nstlim=50000000, 
dt=0.002, 
nrespa=1, 
irest=1, 
nbxe5, 
ntxo=2, 
ntpr=25000, 
ntwr=25000, 
ntwx=25000, 
iwrap=1, 
ioutfm=1, 
ntwprt=0, 
\j1264=1, 
ntt=3, 
temp0=303.15, 
gamma_lIn=1, 
ig=-1, 
ntp=1, 
barostat=2, 
presO=1.01325, 
taup=2, 
ntc=2, 
ntf=2, 
ntb=2, 
cut=9, 
/ 
Se 
6. Building the mutated PAM systems was facilitated by VMD and tleap. In building the NAG system, for example, delete the 2nd G base atoms in the NGG 
PAM first, re-name the original G as A, and then use tleap to complete the missing base atoms of A. In VMD, adjust the orientation of A to avoid potential steric 
clashes if present. The generation of NNGG system followed the same strategy. 
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